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ABSTRACT 

The PRIsm MUti-object Survey (PRIMUS) is a spectroscopic galaxy redshift survey to z ~ 1 
completed with a low-dispersion prism and slitmasks allowing for simultaneous observations of ~ 2, 500 
objects over 0.18 deg 2 . The final PRIMUS catalog includes -130,000 robust redshifts over 9.1deg 2 . 
In this paper, we summarize the PRIMUS observational strategy and present the data reduction 
details used to measure redshifts, redshift precision, and survey completeness. The survey motivation, 
observatio n al tec hniques, fields, target selection, slitmask design, and observations are presented in 
iCoil et aTJ (|201lD . Comparisons to existing higher- resolution spectroscopic measurements show a 
typical precision of er z /(l + z) = 0.005. PRIMUS, both in area and number of redshifts, is the 
largest faint galaxy redshift survey completed to date and is allowing for precise measurements of the 
relationship between AGNs and their hosts, the effects of environment on galaxy evolution, and the 
build up of galactic systems over the latter half of cosmic history. 

Subject headings: galaxies: distances and redshifts, galaxies: evolution, galaxies: high-redshift, large- 
scale structure of universe, surveys 



1. INTRODUCTION 

The PRIsm MUlti-object Survey (PRIMUS) is a new 
spectroscopic survey of faint galaxies focused on mea- 
suring the evolution of galaxy properties and large- 
scale structure since z ~ 1. PRIMUS uses a hybrid 
technique between very low-resolution redshifts mea- 
sured from galaxy photometry and high-resolution spec- 
troscopy. Rather than using a more traditional grating 
or grism, we utilize a specially designed low-dispersion 
prism to observe the full optical spectra of up to — 3, 000 
objects simultaneously. Combined with the 0.2 deg 2 
field of view provided by the IMACS spectrograph on 
the Baade I 6.5m telescope at Las Campanas, PRIMUS 
spectroscopically measures galaxy redshifts in much less 
time than would be required with more traditional high- 
resolution spectroscopic techniques (~ 10, 000 galaxies 
per night) and reaches depths of i — 23 in one hour. 
PRIMUS redshifts have typical precision of a z j (1 + z) = 
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0.005. PRIMUS spectroscopy focuses on fields with ex- 
isting optical imaging for target selection and emphasizes 
areas with existing high-quality multiwavelength imaging 
from GALEX, Spitzer, Chandra, and XMM-Newton. The 
combination of PRIMUS redshifts with this wealth of 
multiwavelength data allows for detailed measurements 
of galaxy colors, luminosities, star formation rates, and 
stellar masses. 

Traditional spectroscopic redshift surveys require con- 
siderable resources to be completed, especially as one 
probes galaxies beyond the local univer se. For exam- 
ple, the DEEP 2 survey (|Davis et al.ll2003|) required more 
than 80 nig hts on Keck and zCOSM OS dLillv eja!1[200l 
and VVDS (|Le Fevre et al.ll2Cl05l IGarilli et al 1120081) uti- 
lized over one hundred VLT nights each. With this large 
amount of observing time, these surveys probe — 20, 000 
galaxies over a few square degrees. PRIMUS represents 
a new paradigm in galaxy redshift surveys which enables 
the efficient survey of large astronomical volumes with 
much less observation time than required with traditional 
higher-resolution spectroscopy. The methods pioneered 
by PRIMUS have already been adopted and refined for 
other pri sm spectroscopic surveys of the high-redshift 
universe ([Kelson et al.ll2012i: [Just et al.|[2012l: iGrav et~aT1 
lin prerJ h 

Science with t he PRIMUS dataset is ongoing. In 
lAird et al.l <l2012Tl . we measure X-ray observations in 
PRIMUS fields to find evidence the presence and ac- 
cretion of AGNs does not depend o n the s tellar masses 
of their host galaxies. iWong et al.l (|2011l) investigated 
tidally triggered star formation in isolated close galaxy 
pairs and measured an enhancement in the specific star 
formation rate due to tidal interactions in these galax- 
ies. We investigated optically "red and dead" galaxies 
which show signatur es of star format ion from their in- 
frared photometry in lZhu et al.l f|201 lh and show that a 
significant fraction of red-sequence galaxies have ongo- 
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ing obscured star formation. In [Moustak as et all (|2012j) . 
we present the stellar mass function of PRIMUS galax- 
ies and conclude that star forming galaxies are quenched 
more strongly with decreasing stellar mass and the ma- 
jority of the stellar mass buildup within quiescent galax- 
ies occurs around ~ 10 10 ' 8 Mq. 

This paper describes the observations, data reduction, 
redshift fitting, redshift precision, and lessons learned 
from our experience with low -resolution prism spec- 
troscopy. A companion paper (|Coil et al.ll20H (here- 
after Paper I) describes the survey design and character- 
istics and provides an overview of the data taken as part 
of the PRIMUS survey. The paper is outlined as follow. 
Section [2] outlines the details of PRIMUS observations 
including details of our nod & shuffle spectroscopy, cali- 
bration data, and prism characteristics. In fj3j we provide 
a detailed discussion of our image reduction and spectral 
extraction including corrections due to scattered light 
and our wavelength calibration. We summarize our mea- 
surements of the flux corrections for our PRIMUS spectra 
in S|4] Section [5] introduces the spectral libraries we uti- 
lize to measure redshifts from PRIMUS spectroscopy and 
detail the fitting algorithm. The redshift precision and 
outlier rates are presented in £ 16.11 and completeness es- 
timates are presented in [J7] before concluding in [J5] The 
first release of PRIMUS redshifts has been completed and 
can be found at http:/ /primus. ucsd.edu. 

2. OBSERVATIONAL STRATEGY 

Traditional wide area multiobject spectrographic sur- 
veys utilize either slit masks or fiber-fed spectrographs 
to obtain spectra for tens or hundreds of galaxies simul- 
taneously. By replacing the traditional grism or grating 
dispersive element with a low dispersion prism, PRIMUS 
utilizes slit masks to obtain spectra for thousands of 
galaxies simultaneously. Figure Q] shows a mosaic for 
all eight IMACS CCDs for one PRIMUS science expo- 
sure illustrating the large multiplexing power available 
with our low-dispersion prism technique. We design each 
PRIMUS slitmask to provide the maximum number of 
high priority galaxies while also ensuring we dedicate 
enough slits on each mask for calibration purposes. 

PRIMUS utilized existing optical photometric catalogs 
for target selection and mask design. Paper I details the 
galaxy and AGN selection criteria, mask design, sparse 
sampling, and photometric catalog details used in our 
mask design. Here, we focus on the observational and 
analytical methods used to obtain spectra for PRIMUS 
targets. 

All PRIM US spectroscopy was observed in n od & 
shuffle mode ()Glazebrook fc Bland-Hawthorn 2001;) with 
IMACS on the Baade Magellan I telescope at Las Cam- 
panas. Figure [2] shows the PRIMUS prism resolution and 
dispersion as a function of wavelength. The prism reso- 
lution is a strong function of wavelength and has lowest 
resolution in the red. Robust sky subtraction is vital for 
PRIMUS redshift fitting as most galaxies are fit based on 
the shape of the continuum (galaxies with strong emis- 
sion lines include those lines in the redshift fitting, but 
many galaxies in our sample have weak or absent emis- 
sion lines). As nod & shuffle allows us to measure the 
sky flux in the same pixels as we measure object flux, we 
can perform robust sky subtraction. 

Rather than nodding the mask such that objects moved 
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Fig. 1 .— Mosaic of all eight IMACS CCDs for a single PRIMUS 
science exposure. Each trace of light corresponds to two individ- 
ual slits which have been shuffled on the detector, resulting in 4 
traces per object. The very bright traces in the image correspond 
to alignment stars utilized to ensure that the slitmask has been 
centered properly when pointing the telescope. 
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Fig. 2. — Resolution and dispersion of the PRIMUS prism versus 
wavelength. The resolution and dispersion generated by the low 
dispersion prism is a strong function of the wavelength with in- 
creased resolution in the blue but low resolution at 9000A — 1/im. 
The low resolution allows us to observe the full optical spectrum 
of a target galaxy on only ~150 pixels on the detector and up to 
~3000 galaxies on a single mask. 
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along one long slit, we adopted a strategy of cutting two 
slits (often referred to as slit A and slit B throughout 
this paper) for each object on the mask separated by 
(Aa, AS) = (-1'.'6, -3'.'2) when we designed 0'.'8 wide 
slits and (Aa, AS) = (-2'.'0, -4'/0) when we designed l'.'O 
slits. By nodding in both directions on the sky rather 
than in a single direction, we minimize the impact of a 
bad column on the final extracted spectrum of an ob- 
served galaxy; a galaxy moved in both directions will 
not lose all of the information from the wavelength asso- 
ciated with the bad column. During each telescope nod, 
the detector was also shuffled 8 (10) pixels for masks 
with 0'.'8 (l'.'O) wide slits. An example of nod & shuffle 
traces for a single object is show in Figure [3] Due to nod- 
ding between two slits, the final observation includes four 
traces for each object; two traces with only sky illumi- 
nation and two with both sky and object contributions. 
In November 2007, modifications to the IMACS instru- 
ment lead to a slight change in this layout; objects after 
this date were observed in the outer two traces and the 
inner two traces contained only sky light. Throughout 
this paper we will often discuss slits (either object slits 
or sky slits), referring to the spectrum observed through 
a single slit (in either object or sky position). When we 
discuss the "full object trace" we refer to the collection of 
the 4 traces shown in Figure |3] including the 2 sky traces 
and 2 objects traces coming from slit A and slit B in nod 
& shuffle mode. 

Science frames were typically observed in four 16 
minute nod & shuffle exposures for a total exposure time 
of 64 minutes (32 minutes in each nod & shuffle object 
slit). In COSMOS, we exposed deeper due to the wealth 
of high-quality data in the field; COSMOS masks in- 
cluded six 16 minute exposures for a total of 96 minutes 
per mask. 

3. IMAGE ANALYSIS AND SPECTRAL 
EXTRACTION 

We utilize a multi-stage process to first remove contam- 
ination light from the the full two-dimensional science 
frames and then to extract the one-dimensional spec- 
trum of each observed object. In this section, we provide 
the details relevant to our processing of IMACS images 
and the methods we utilize to obtain high-fidelity low- 
dispersion spectra via nod & shuffle sky subtraction. Fi- 
nally, we discuss the modeling methods used to extract 
object spectra. 

3.1. Large Scale Scattered Light 

IMACS images suffer from large diffuse halos of light 
around bright sources. The total light contained in this 
halo is correlated with the flux of nearby objects and is 
extended enough to influence the light gathered in nearby 
object traces. The scattered light crosses chip gaps, in- 
dicating that the halos are not due to scattering within 
the CCDs themselves but are likely created in the cam- 
era optics. The presence of excess light will lead to errors 
in sky subtraction and can easily contaminate extracted 
spectra if not carefully corrected. 

Figure H] illustrates the presence of the large scale halo 
light near two object trace sets from a PRIMUS science 
exposure. On the right, we show two slices through the 
image; the red and blue slices are taken at the locations 
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Fig. 3. — Example of a full bright-object trace from PRIMUS 
spectroscopy. PRIMUS observations are completed in nod & shuf- 
fle mode with the object nodded between two slits separated in 
both right ascension and declination. This pattern results in four 
traces on the final image, two traces of pure sky light and two 
traces of object combined with sky. In this figure, red is to the left 
and blue to the right. For reference, the bright sky line near pixel 
73 (83) in the bottom (top) trace corresponds to 5577 A. The res- 
olution and dispersion of the PRIMUS prism are strong functions 
of wavelength; the forest of emission lines in the red portion of the 
spectrum are blended at PRIMUS resolution. 

of the vertical red and blues lines in the two-dimensional 
image. The clear signal of light between the trace sets is 
the signature of the large scale scattering halo we char- 
acterize and remove using the details in this subsection. 

In order to remove the contribution from the large- 
scale scattering light component, we first stack two- 
dimensional images of full traces with no targeted object 
(originally milled to characterize sky emission) and use 
this stacked image to create a model for the scattered 
light profile. We parameterize the scattered light term 
as a convolution of our image with kernel, k(x, y) 

k(x,y)= (l+^+J^) " (i) 

where x and y correspond to the pixel coordinates on 
the detector and r and p are free parameters we fit to 
characterize the shape of the scattered light halo. 

We create a two-dimensional model for the scattered 
light by coadding all of the science exposures for a given 
pointing to increase the signal-to-noise ratio in the scat- 
tered light and remove cosmic rays. We mask the lo- 
cation of full traces on the coadded image and then 
fit for the kernel scale length, ro, and power-law in- 
dex, p, over the full 0.2 deg 2 IMACS focal plane using 
the Le venberg-Marquard t method with the IDL routine 
mpf it (|Markwardtll2009l ). 

We find that a simple convolution of the best-fit kernel 
and the original science image is not sufficient to pro- 
vide a full correction to the scattered light due to several 
factors. First, the amplitude of the scattered light with 
respect to the incident light from the slits varies as a func- 
tion of position on the focal plane. Second, we find that, 
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Fig. 4. — Left: Two dimensional sub-region from one PRIMUS 
science exposure. Two object trace sets (each object set is com- 
posed of four traces) are shown. The presence of light between the 
two trace sets is the signature of the large scale scattered light halo 
we remove before performing object extraction from our data. The 
red and blue vertical lines illustrate the locations where two slices 
were made through the image. The slices show that this scattered 
light contributes several hundred counts per pixel and is a clear 
contaminant to our spectra. Furthermore, since the scattered light 
has two dimensional structure on our images, performing nod & 
shuffle sky subtraction without removing the contaminating signal 
will result in poor sky subtraction. 

in detail, the large-scale halo is not circularly symmetric 
as assumed by equation (JT|). In order to obtain a better 
model of the scattered light, we define two non-circular ly 
symmetric kernels based on the best fitting k(x, y) 

9 = t&n- 1 {y/x) (2) 
h(x,y) = k(x,y)cos(6) (3) 
k 2 (x,y) = k(x,y)cos(6 + iT/2). (4) 

We solve for the spatial dependence of the halo term 
by convolving the masked image with the three kernels 
to generate three halo terms 

A/halo = I ® k (5) 

-Dl.halo =I®ki (6) 
£>2,halo = I®k 2 . (7) 

We then solve for the spatial variation in the halo by 
creating the final halo image, /halo, and solving for Wi to 
minimize the light between object traces 

Jhaiofo y) = (wo + wix + w 2 y + w 3 x 2 + Wiy 2 ) + (8) 
Afhalo(«>5 + w 6% + w 7 y + w$x 2 + w 9 y 2 ) + 
u>io-Di,halo + wii-D 2 , halo- 
When applying the best fitting halo image, /halo, con- 
structed in this manner, we normalize the counts in the 
halo to the median counts in each individual exposure 
to adjust for exposure to exposure differences. Figure 
[5] follows Figure |4] but illustrates the same two dimen- 
sional image after the large scale halo model has been 
subtracted. The slices shown on the right in the figure 
highlights the success of this technique; little light re- 
mains between the object traces. 



3.1.1. Small Scale Scattered Light 

In addition to the large scale scattering halo, our 
IMACS spectroscopic images have a second small-scale 
scattering component on 5-10 pixel scales. This halo 
component is most apparent in the red portion of the 
spectra and is non-axisymmetric with respect to the light 
in the full traces. On average, this small scale halo light 
corresponds to approximately 1% of the flux of the sky 
While the cause of this excess light is unknown, the local 
nature of the light suggests it may be due to diffusion of 
red light in the CCDs. 

For each CCD, we first locate all full trace sets with- 
out another trace within 20 pixels. For each night of ob- 
servations, we stack the two dimensional image of these 
isolated traces after performing a nod & shuffle sky sub- 
traction. We then quantify the presence of excess light 
outside of the sky traces in the composite image. We 
assume the object trace experiences the same amount of 
scattering as the sky trace (although mirrored across the 
trace) and construct a final empirical model for the small 
scale scattered light in the spectra. In the top panel of 
Figure we show an example of such a stacked image. 
The scattered light is apparent especially in the red (left 
side). We take the scattering light in the 7 (9) pixels (for 
0.8" (1".0) wide slits) above the top trace and that in the 
7 (9) pixels below the bottom trace, and reproduce the 
light for each trace based on the same nod & shuffle mode 
as observations. In the bottom panel of Figure|51 we show 
the empirical model constructed from the stacked image 
in the top panel. Finally, we duplicate this model for 
each trace on the CCD to subtract the scattering light 
from the image. Overall, this correction is ~ 1% of the 
sky and thus an important contamination to remove from 
our data before extracting the one-dimensional spectra. 

3.2. Flatfield Correction 

At this point, we have performed no flat fielding of 
the PRIMUS two-dimensional images. When fitting 
PRIMUS spectra to obtain redshifts, we are not con- 
cerned with the normalization of the spectra (which is 
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Fig. 5. — Two dimensional subregion shown in Figure[4] after the 
best fitting model for the large scale halo has been subtracted from 
the image. The slices through the image plotted on the right do 
not show signal between the traces. 
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affected by vignetting and slit losses) but removal of 
pixel-to-pixel changes from the detector is vital in order 
to not introduce false spectral features in our extracted 
spectra. When constructing a flat field correction for our 
data, we first remove any large-scale gradients present in 
our observed twilight sky flats. The resulting flat field 
measures the pixel-to-pixel variations across the detec- 
tor. Since the twilight flat is observed as a single image 
and our science frames are observed with nod & shuffle, 
we cannot simply divide the science data by this pixel- 
to-pixel flat field correction. 

Using our normalized flat-field, we construct a "nod 
& shuffle" corrected flat field by finding the mean be- 
tween the observed flat field image and one shifted by 
our shuffle length. In general, PRIMUS spectra are sky- 
dominated, so the assumption that each nod position 
contributes evenly to the total contribution to the flat 
field is sufficient to correct all but the most severe blem- 
ishes in the flat field. Areas on the detector that are 
affected by flat-field corrections of more than 25% are 
masked as possible bad pixels and are ignored through- 
out the redshift fitting process to ensure that these large 
flat-field corrections, where our simple assumption may 
not be ideal, do not affect our final spectra. Less than 1% 
of the pixels are masked to this large flat fielding error. 
This final flat frame is utilized in our forward modeling 
of the PRIMUS science frames; in contrast to more tra- 
ditionally approaches, we do not directly apply the flat 
field correction to our images. 

3.3. Spectral Modeling & Extraction 

Having corrected for the scattered light in the two 
dimensional images, we are ready to extract the one- 
dimensional object spectra from the images. In 
PRIMUS, we utilize a forward modeling approach to 
spectral extraction; we first perform an extraction with- 
out information on the spatial profile of each object. We 
then perform more robust extractions to construct a one 
dimensional profile and perform an extraction based on 
that profile to obtain the final extracted spectra. This 
section details the steps involved with calibrating the two 
dimensional images and extract science-quality spectra. 

3.3.1. First pass extraction 

We perform spectral extraction of the halo-corrected 
two dimensional PRIMUS spectra in a three pass 
method. The first pass of extraction allows us to derive 
the spatial profile used to extract each object, the second 
pass refines the sky subtraction and the final pass mea- 
sures the one-dimensional extracted spectrum for each 
object. 

In the first pass, we use a simple model for the sky and 
object light in order to build a more detailed spatial pro- 
file for each object in the subsequent extraction passes. 
Extraction is performed in a forward-modeling sense on 
a column-by-column basis in the images. For each col- 
umn, we identify pixels associated with each PRIMUS 
trace. We construct models, Si, the sky light in slit i 
(cither slit A or slit B), and Oi, the object light in cither 
slit position. In the first pass extraction, we construct 
simple models for each 

JVpix 

Si = "'-..A-// ( 9 ) 
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Fig. 6. — Top: Stacked image of all of the isolated traces on a 
single PRIMUS observation which illustrates the presence of small- 
scale scattered light (for examples the "spurs" of light extending 
above and below the traces). Bottom: Empirically constructed 
model used to r emove the small-scale scattered light constructed 
as described in t|3.1.1l 

AW 

O i = Y i w jS j f j (10) 

3 

where fj is the flat field contribution to pixel j in the 
trace, Sj are delta-functions in each of the pixels across 
the trace, and w s j and w a j are the fluxes measured in 
sky and object pixels prior to the flat-field contribution. 
Before fitting, we construct a model for the full four trace 
contribution to the column from the object; the sky con- 
tributes flux only in the sky slits and both the sky and 
object model contribute flux in the object traces. Finally, 
this model is convolved with a boxcar with the width of 
the mask slit and a Gaussian with a 1.5 pixel width to 
represent the instrumental profile. By solving for w s .j 
and w j, we construct a rough model for the sky and 
object illumination in the slit. 

3.3.2. Second and third pass extractions 

The second pass of extraction uses the output from the 
first pass in order to more robustly extract the object 
spectra. We start this process by generating a spatial 
profile of the object light in each slit. We sum the object 
light, as measured in the first pass, along the spectral 
direction of the PRIMUS trace to measure pj , the total 
light in the slit as a function of pixel j across the slit. We 
do this sum separately for each A and B nod & shuffle 
position; in an ideal observation, these spatial profiles 
would be identical, but due to guiding errors or errors in 
the mask cutting, the objects may not be centered at the 
same pixel in both positions. 

We fit pj with a Gaussian model (solving for the cen- 
troid and Gaussian width) in order to model the spatial 
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Fig. 7. — Results of the second pass of PRIMUS spectral ex- 
traction. The black data points with error bars shows the actual 
PRIMUS data while the red (dot-dashed) and grey denote the best 
fitting object and sky contributions to the four traces. The final 
combined model is shown by the grey line. 



profile of the object. While this measurement is robust 
for well-detected objects in a single exposure, the low- 
est signal-to-noise objects have derived profiles that are 
heavily influenced by noise in the images. Rather than 
using noisy profiles for the low signal-to-noise objects, we 
force these objects to have the median width and centroid 
of the other objects on the same CCD. Based on the in- 
ferred centroid, fii, and profile width, Cj, of object i, we 
construct a new object model 
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We perform the extraction column-by-column using 
this spatial profile for the object traces rather than the 
S- function formulation in equation 1101 Figure [7] shows a 
slice through one relatively bright PRIMUS object trace 
set. The red lines show the best fitting object profile 
while the blue shows the contribution from the sky. The 
grey line shows the sum of the sky and object models 
compared to the black data points with errors showing 
the actual PRIMUS data. We utilize the best fitting sky 
parameters from this pass to create a two-dimensional 
sky image which is subtracted from the halo-corrected 
two-dimensional science frame. We perform a final ex- 
traction utilizing the spatial profile derived above to ex- 
tract e ach object tr ace with an optimal extraction algo- 
rithm (|Hornd 11986) which maximizes the signal-to-noise 
ratio in the extracted spectra compared to a simple box- 
car extraction. 

3.4. Wavelength calibration 

We perform wavelength calibration using a three step 
process. Initially, we apply a model wavelength vec- 
tor to each extracted object and sky spectrum based 
on the optical model for the IMACS spectrograph and 
PRIMUS prism. In this model, the pixels are spaced 
roughly evenly in A~ 3 . Figure [8] shows the wavelength as 




Fig. 8. — Wavelength as a function of pixel position along the 
trace for PRIMUS observations. This characteristic curve has pix- 
els roughly evenly spaced in A -3 . When solving for the final wave- 
length vector for each object in the survey, we measure the zero- 
point offset between this fiducial model and the observations based 
on the locations of key sky emission features observed in each sky 
slit in our nod & shuffle exposures. 



a function of pixel from this optical model. In order to 
ensure that this model correctly constrains the final per- 
formance of the prism, we observed images illuminated 
by helium arc lamps on the flat field screen. We per- 
form a simple extraction of the arc images and compare 
the spacing of the well-separated helium lines with the 
expected positions based on the prism dispersion model 
and find excellent agreement for all PRIMUS observing 
runs. Figure |9] shows one such comparison; throughout 
all runs, the general form of the relationship between 
pixel and wavelength remained fixed. 

While the optical model of the instrument and prism 
perform well in predicting the dispersion in PRIMUS 
spectra, it cannot accurately predict the zeropoint for 
the wavelength solution on the scales needed to achieve- 
ment robust redshifts. While the PRIMUS spectra are 
low resolution, the brightest sky lines are still readily 
detected and make ideal zeropoint normalizations to our 
wavelength solutions. Based on the strong night sky lines 
at 5577.34, 5895.00, 6300.3 A, we derive a zeropoint offset 
for each slit on a CCD. Due to the low resolution nature 
of our spectra, we find the most robust zeropoints are 
determined by doing a cross correlation between a model 
sky spectrum and our observed sky in windows around 
each of the bright sky lines. Based on this cross corre- 
lation, we measure an offset for each object on a given 
PRIMUS exposure. In order to properly calibrate ob- 
jects where a portion of the bright sky lines are affected 
by bad columns and to protect against contamination 
from cosmic rays which may occur if we apply a slit-by- 
slit zeropoint correction, we fit a low order polynomial to 
the the final zeropoint shifts as a function of position on 
the IMACS focal plane. We then apply the correspond- 
ing model value to each slit in our science observation. 

The zeropoint derived from the sky lines is appropri- 
ate if the slit is filled evenly by the object or for slits 
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Fig. 9. — Comparison of the wavelength solution determined by 
using an optical model for the IMACS optical path combined with 
our refinement based on the location of prominent features in the 
sky spectrum with an extraction of a helium lamp observation 
through the same slit. Using a technique based on sky line features 
allows us to determine wavelength corrections based on simultane- 
ous observations with our object spectra rather than extrapolating 
corrections from images which bracket our science frames. The lo- 
cation of key lines in the helium spectrum are highlighted in by 
dot-dashed lines. 

in which the object is centered. Due to alignment er- 
rors, mask cutting errors, astrometry errors in the pho- 
tometric catalog, etc, this assumption is not always true 
for PRIMUS observations. Ignoring the effect of non- 
aligned slits can lead to systematic errors in our final 
wavelength calibration, and thus we utilize the ~ 30 F- 
stars observed on each PRIMUS mask in order to derive 
systematic zeropoint offsets in the object spectra after 
the sky-line normalization has been ap plied. For each 
target F-star candidate, we fit a grid of iKuruc^ (|1993l ) 
stellar models to the stellar photometry in order to deter- 
mine an appropriate model. We convolve the best fitting 
stellar spectrum and measure, through a x 2 search, any 
residual offset between the model and stellar wavelength 
model. Typically, the final shifts are small (< 0.1 pixels) 
and vary smoothly as a function of position on the focal 
plane (as one would expect if the shifts were due to align- 
ment errors or large scale drifts in the astrometry). We 
fit the shifts as a function of position on the Baade focal 
plane and apply the zeropoint correction to each object 
spectrum. 

3.5. Coadding Extracted Spectra 

Throughout the extraction and wavelength calibration 
process, we work on each individual science exposure sep- 
arately rather than on a coadded image of all of our ex- 
posures. Primarily, we do this to protect against guiding 
errors, which can accumulate from exposure to exposure, 
leading to systematic problems when trying to extract a 
coadded image. Furthermore, as we work in nod & shuf- 
fle mode, if a subset of the exposures were subject to 
changes in the sky brightness on timescales shorter than 
our nod time, these images will have systematically poor 
sky subtraction; including them before extracting will 



lead to systematic errors in our extracted spectra. 

Before coadding individual exposures, we first search 
for frames which do not have the quality desired to be 
included in the final stacked science data. If a frame was 
observed under conditions in which the nod & shuffle pro- 
cedure fails due to sky conditions changing on timescales 
shorter than our nod timescale, the sky subtraction will 
leave residual light (or deficient light) in one object slit 
compared to the second. This will manifest itself both 
in systematic errors when the difference between the ob- 
ject spectra in slit A are compared to slit B, Oa — Ob- 
Furthermore, remnant sky light in the extracted spectra 
will lead to systematically different signal-to-noise ratios 
in slit A, SNa, versus that in slit B, SNb- We isolate 
any exposures with significant difference from the mean 
in the < Oa — Ob > versus SNa/SNb plane and iden- 
tify those exposures as possibly contaminated exposures 
and reject those exposures when constructed the coad- 
ded final spectra. Of the 705 science images we observed 
for PRIMUS, 12 exposures were rejected based on these 
criteria. 

After any problematic exposures have been removed, 
we coadd the extracted one-dimensional spectra mea- 
sured in each of our exposures to determine the final spec- 
trum for each object. In order to weight by the quality 
of the data in each exposure (for example due to varying 
airmass or transparency between exposures), we weight 
each exposure according to the mean signal-to-noise ra- 
tio of the mask. We construct this mean signal-to-noise 
ratio by fitting the signal-to-noise of each object spec- 
trum as a function of object magnitude and determine 
the value of this trend at i = 21. We do not coadd the 
extracted spectra from slit A and slit B in each expo- 
sure. If we were to coadd these spectra, the fact that 
our objects are shifted both vertically and horizontally 
on the IMACS focal plane means we would need to shift 
the two measurements spectrally in order to coadd and 
would result in the noise in our final coadded spectrum 
being correlated. While this correlation could be tracked 
and treated properly in the redshift fitting, we preferred 
to not coadd the two slits and simultaneously fit the spec- 
trum from slit A and slit B when solving for redshifts. 

After our full extraction, scattered light corrections, 
and nod & shuffle sky subtraction, we find that slits 
drilled without any objects in them show small-scale 
residual light. In order to correct for this residual 
sky-subtraction error, we coadd the extracted counts in 
empty slits separately in slit A and slit B for each of 
the eight IMACS ccds for an entire night of observa- 
tions. This results in 16 wavelength-dependent correc- 
tions for each night (one for each of slit A and slit B 
on the 8 CCDs). We interpolate this mean residual to 
the wavelengths of each object in our coadded spectra 
and subtract it before proceeding with redshift fitting. 
Figure [TO] shows the magnitude of this effect compared 
to scaled sky spectra. Overall, this effect is 1-2% of the 
sky flux and thus an important correction. Removing 
this term greatly improves the final agreement between 
the extracted counts in slit A compared to slit B and 
redshift failures rates for galaxies with known redshifts. 
On nights in which, primarily due to hardware problems 
with the mask cutting, we had too few empty slits to 
make these measurements we utilized the derived correc- 
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Fig. 10. — Magnitude of the residual light correction, derived by 
coadding the extracted spectra from slits in which no objects were 
observed, as a function of wavelength. The orange (dashed) and 
blue (dot-dashed) lines show the median correction for slit 1 and 
slit 2 respectively. The small scale correction contains about 1-2% 
of the flux of the sky and thus is an important correction when 
studying faint objects. 

tion from the nearest night with viable data. 

4. SPECTROPHOTOMETRIC CALIBRATION 

Accurate spectrophotometry is essential for determin- 
ing reliable redshifts for PRIMUS because our spectral 
resolution is too low to rely on resolve individual nar- 
row emission lines; therefore, for many objects the red- 
shift is largely determined by the spectral shape (e.g., the 
strength of the 4000- A break). We emphasize that our 
objective is to obtain reasonably accurate relative spec- 
trophotometry; the absolute normalization of each spec- 
trum, which depends on the amount of light lost from 
the slit due to the physical extent of each object relative 
to the slit width, variations in the PSF, pointing errors, 
and so forth, is not relevant for determining the redshift 
with our method. 

Traditionally, standard-star observations are used to 
derive a sensitivity function to flux-calibrate the observed 
spectra. By contrast, we approach the problem of spec- 
trophotometric calibration from a forward modeling per- 
spective. In other words, we apply our best estimate of 
the wavelength-dependent throughput to the model tem- 
plates used for redshift fitting (see $5]). Thus our prin- 
cipal goal is to determine the throughput as a function 
of both position on each CCD and time; we then employ 
the throughput estimate when fitting for redshifts. 

We first utilized slitless DA white dwarf observations 
taken in Jan 2006 to obtain an overall throughput curve 
for the IMACS instrument and PRIMUS prism as shown 
in Figure [TT] We are able to calibrate the wavelength 
of the extracted slitless spectrum from the broad B aimer 
lines present in the DA spectrum. Next, we correct the 
spectra for atmospheric extinction at the mean airmass 
of each mask utilizing the Gemini South atmospheric ex- 
tinction curve. Clearly one correction curve obtained at 
the beginning of the survey is not appropriate to fully 
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Fig. 11. — Baseline IMACS throughput used as a zeroth-ordcr 
flux calibration term throughout the PRIMUS flux calibration. 
This curve was obtained by observing a DA white dwarf on CCD 
5 of the IMACS array. The broad absorption line features allowed 
us to perform wavelength calibration on the slitless spectrum. We 
utilize the observed sky emission in each PRIMUS slit to extend 
the flux solution between each of the IMACS CCDs and tie each 
CCD to the same flux solution. We finally apply a correction based 
on the ~ 30 F-stars we observe on each mask to remove any varia- 
tions from the overall IMACS response function over the course of 
the survey's operations. 

calibrate our data in subsequent runs; we utilize this nor- 
malization as a baseline. More detailed corrections are 
applied to the data based on observations of F-stars and 
the night sky to remove the higher order fluxing residuals 
from this baseline correction. 

In order to quantify how the flux correction varies be- 
tween each of our PRIMUS observing runs, we utilize 
two tracers for the sensitivity measured on each of our 
PRIMUS masks. We utilize the ~ 30 F-stars on each 
mask to measure an average overall sensitivity function 
for each PRIMUS run and use the sky spectra in each 
mask to correct for any variations in the inter-ccd sen- 
sitivity corrections. While, in practice, we could gen- 
erate mask by mask corrections, the variations between 
masks in a given run are small; averaged corrections al- 
low the use of more F-star and sky measurements when 
constructing our fluxing terms leading to higher signal- 
to- noise ratio corrections. 

The baseline IMACS sensitivity is tied to CCD 5 in the 
IMACS array. In order to normalize each of the other 
CCDs, we utilize the sky spectrum observed in the "off" 
position in our nod & shuffle observations. Based on the 
assumption that the sky spectrum is constant over the 
IMACS field of view, we measure the mean sky spectrum 
on each of the eight IMACS CCDs. We then measure the 
ratio between the mean sky measured on CCD 5 to the 
mean sky measured on each of the other CCDs. Figure 
[T"2l shows this procedure for one PRIMUS mask. The 
black spectra show the CCD 5 mean sky spectrum and 
the blue shows the mean for each of the other CCDs. The 
red line is a bspline fit utilized to remove the variation 
between IMACS ccds. Note that we assume the sky spec- 
trum does not vary over the IMACS field of view. While 
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this is not an unreasonable assumption for the sky con- 
tinuum, emission line strengths can vary on arcminute 
scales. The grey lines in Figure [12] illustrates this issue. 
As we ultimately fit a smooth bspline to the data, the 
few percent variations in the sky emission lines do not 
significantly affect our methodology. We do not expect 
the relative normalizations between the CCDs to be a 
strong function of time; by combining data over multiple 
masks and nights allows us to construct the highest possi- 
ble signal-to-noise ratio correction. After examining the 
time dependence of these corrections, we chose to create 
these "inter-ccd" corrections for each month of PRIMUS 
observations. Figure [TBI shows the variation with time of 
these correction terms. Overall, the curves for each CCD 
follow the same general shapes over the full course of the 
PRIMUS observations. Some differences arise likely due 
to modifications to the instrument and electronics over 
the course of our observations. 
With each of the CCDs on the same spectrophotomet- 
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Fig. 12. — Derivation of the inter-ccd spectrophotometric cor- 
rections for one PRIMUS mask. The baseline IMACS through- 
put and comparison to F-stars observed on PRIMUS masks ties 
the PRIMUS spectra spectrophotometrically to IMACS CCD 5. 
Based on observations of the sky spectrum while in the "off" posi- 
tion in our nod & shuffle observations, we correct for the response 
differences between CCDs by constructing the mean sky spectrum 
on each CCD (blue). The ratio between this mean spectrum and 
the mean sky spectrum from slits on CCD 5 (black) gives a spec- 
trophotometric correction to bring each ccd to the same system 
(grey). Lastly, we fit the observed ratio with a smooth bspline to 
characterize the low-order response variations between CCDs (red) . 
As the temporal variation of these corrections is low, we combine 
observations over each month of PRIMUS observations when cal- 
culating the final inter-ccd flux correction to provide the highest 
signal-to-noise ratio correction possible. 
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Fig. 13. — Variation of the inter-ccd flu x corrections described 
in the text and illustrated in Figure 1121 as a function of time. 
For each CCD, we show the final mean correction used to correct 
for variations between CCDs for each month of PRIMUS observa- 
tions. Overall, the general form of the correction remained constant 
throughout the survey, but modest changes occur primarily due to 
instrument maintenance throughout the course of the survey. 

ric system, we next bring the total spectrophotometric 
system of each mask onto the final survey system. While 
the baseline IMACS correction corrects for a general re- 
sponse function of the instrument, this response func- 
tion can vary over time. In order to remove this vari- 
ation, we utilize the ~ 30 F-stars we observed on each 
mask as our spectrophotometric standards. For each F- 
star, we fit its broadband photometry against a model of 
Kurucz (1993) stellar atmosphere models. We measure 
the ratio between the observed PRIMUS spectrum and 
the best-fitting stellar model which has been convolved 
to PRIMUS resolution. Figure [JJ] (a) shows this fitting 
process; the datapoints with errorbars sho w the PRIMUS 
data and the grey line is the convolved iKuruczl (|1993[ ) 
model. The resulting correction vector is the ratio be- 
tween the observed and model data. For each PRIMUS 
observing season, we combine each of these correction 
vectors to derive a median correction vector for the en- 
tire run as shown in Figure [TH (b) . After careful com- 
parisons, we found that the fluxing corrections within 
a single observing season were consistent, but correc- 
tions between seasons shows considerable variation. Fig- 
ure Q~J] (c) shows the resulting flux corrections for each 
of the PRIMUS observing seasons. The changes illus- 
trated in the figure are often physically motivated. Sev- 
eral major modifications were completed to the IMACS 
instrument throughout the course of our observing pro- 
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Fig. 14. — Derivation of overall flux calibration with PRIMUS F- 
stars. Panel a: The ratio of a PRIMUS-observed F-star (black dat- 
apoints) and model F-star spectrum convolved to PRIMUS resolu- 
tion (grey) is used to create the correction needed to bring PRIMUS 
spectra to a standard system. Panel b: To obtain the highest 
signal-to-noise ratio measurement of the flux calibration, we coadd 
the individual corrections over PRIMUS observing seasons to cre- 
ate the final mean correction. The raw dat a show the me asured 
ratio between observed spectrum and model Kurucz ( 1993) spec- 
trum. The final mean correction (red dashed line) is obtained by 
finding the mean of the individual corrections. The blue dot-dot- 
dot-dashed line shows the ltj region around the mean. Panel c: The 
temporal variation of the F-star flux calibration vectors. Through- 
out the course of the survey, the IMACS team made several im- 
provements and performed regular maintenance to the instrument 
- the results are clear in the increase blue sensitivity of the IMACS 
CCDs throughout the course of the survey. 



gram to increase the throughput and sensitivity of the 
instrument. The corrections clearly shows that these ad- 
justments made sizable improvements, especially at blue 
wavelengths. 

All of these corrections correct the spectral shape of the 
PRIMUS spectra. The main source of small-time varia- 
tion in the flux correction of our spectra which is not re- 
moved in this analysis is due to transparency variations 
throughout each night. These transparency variations 
will drive the overall normalization of our spectra high 
or low but will not affect the spectral shape of the spec- 
tra. When fitting our redshifts, as described in Hb.'Sl we 
tie the normalization of our final model to the broadband 
photometry and allow the normalization of the spectra 
to be set by the photometric data, so transparency vari- 
ations do not affect the quality of our redshift fits. 

5. DETERMINING REDSHIFTS 



5.1. Galaxy Spectral Templates 

When fitting PRIMUS spectra, we chose to compare 
each PRIMUS spectrum to an empirical library of galaxy 
spectra in order to measure the best fitting galaxy red- 
shift for each object. We use empirical templates, rather 
than a more general non-negative linear combination of 
spectral models. At PRIMUS resolution, fine spectral 
features cannot be utilized to help break degeneracies 
between two models which may produce similar low- 
resolution spectra but significantly different when viewed 
at higher resolution. Thus, relying on empirical galaxy 
templates allows us to fit our PRIMUS sample while pre- 
venting the best fitting models from occupying areas of 
parameters space which have no known counterparts in 
reality. 

In order to build our empirical library of galaxy tem- 
plates, we use the AGN and Galaxy Evolutio n Survey 
(AGES) galaxy sample (IKochanek et alJlMtl . AGES 
observed 18163 galaxies with a median redshift of z ~ 0.3 
in the NOAO Deep Wide-Field Survey Bootes field. This 
sample is an ideal source for our empirical templates due 
to the high quality spectroscopy to I = 20.4 and at higher 
redshifts than probed by lower-redshift surveys such as 
SDSS. 

When constructing the sample of AGES galaxies to be 
considered in our library building, we limit the sample to 
galaxies with 0.05 < z < 0.35 and only keep spectra with 
a mean signal-to-noise ratio of at least 4 per pixel in the 
AGES spectrum. This redshift range ensures that Ha fell 
in the observed AGES wavelength range. Furthermore, 
to ensure high-quality spectrophotometry and to create 
a statistically complete sample, we only consider objects 
with 14.5 < /yoga < 19.5. We also only included AGES 
observations that had robust spectrophotometry. These 
quality cuts leave a sample of 3244 galaxies which will 
refer to as the "AGES parent" sample. 

As the observed spectra from AGES do not extend suf- 
ficiently into the blue to cover the full observed PRIMUS 
wavelength range when shifted to higher redshifts, we 
instead use a best-fitting spectral model for the parent 
sample galaxies when deriving the empirical library. For 
each galaxy in the AGES parent sample, we perform a 
detailed spectral modeling of the AGES spectra simulta- 
neously with the measured optical, GALEX, and Spitzer 
photometry. We utilize th e continuum and e missio n line 
fits to AGES galaxies from lMoustakas et al.l (|2011| ) when 
constructing the PRIMUS galaxy template set; the full 
details of the methods inv olved in the fitting of the AGES 
spectra can be found in iMoustakas et alj ([2011D . The 
best fitting parameters were then used to construct a 
full UV to IR model for each AGES galaxy which is used 
throughout the rest of this analysis. 

In order to construct a library of empirical templates 
from the best fitting models, we create a subsample ("ba- 
sis set" hereafter) of the full AGES sample which, when 
convolved to the PRIMUS resolution, fully span the par- 
ent sample. To construct this sample, we first convolve 
each AGES sample galaxy to PRIMUS resolution assum- 
ing it was observed at z = 0.25. We next compute a 
3244x3244 matrix from these convolved spectra; each el- 
ement, Cij, in this matrix is the x 2 value calculated as- 
suming galaxy i could be described by a simple scale 
factor multiplication of galaxy j in the AGES sample. 
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One pair of spectra can represent each other in our basis 
set approach if the \ 2 is smaller than a given maximum 
tolerance Xmax- Based on this matrix, we strive to de- 
termine the minimal number of AGES sample galaxies 
that will span the entire AGES population. 

When constructing the PRIMUS basis set, we set 
Xmax = 2- We then set each element in a new 3244x3244 
matrix to 1 if xfj < Xmax an d if the calculated x 2 was 
larger than our threshold value. We then minimize the 
number of columns subject to the condition that each 
row must have at least one element larger than 0. Math- 
ematically, we solve 

min^Xj (12) 

3 

subject to 

^otijXj>\ (13) 

3 

where Xj — 1 if column j is a basis galaxy, Xj = other- 
wise, and (Xij is the element in our binary matrix. This 
method results in 206 basis set galaxies. 

While the galaxy models utilized to construct this ba- 
sis set is complete in galaxy properties over the range of 
properties probed by the parent sample of AGES galax- 
ies, the redshift limitation when constructing the AGES 
parent removes any information about the 2 700 A Mgll 
emission line strengths from our basis set. To remedy 
this, we isolated any basis set galaxies with 01 (w — g) < 
1.0 and EW[oii] < 15 to be possible Mgll emitters based 
on AGES galaxies that have detected Mgll emission at 
0.4 < z < 0.9 and added a second basis galaxy to the set 
identical to the first with the exception of the addition 
of a 25 A equivalent width Mgll line. This duplication 
is valid as we are simply seeking a set of galaxies which 
span the properties of galaxies at PRIMUS resolution 
but we do not require the minimal set of such galaxies. 
Figure [15] shows 15 randomly selected galaxies from the 
basis set used when fitting PRIMUS galaxies. As the 
largest source of spectral diversity in the galaxies, once 
convolved to PRIMUS resolution, is due to the relative 
strengths of optical emission lines, the basis set is pri- 
marily composed of blue star forming objects, while the 
red galaxy population is well sampled with only a small 
number of basis objects. 

While our final basis set spans the full range of galaxy 
types in AGES when observed at PRIMUS resolution, 
some galaxies in this basis are much more common than 
others. In order to quantify this for later use in the 
PRIMUS redshift fitting, we utilize the basis set to fit 
the full AGES parent sample. For each basis galaxy, we 
count the number of parent galaxies that were best fit 
by it and record this value as the basis weight of the 
template, 6j. Figure \TE\ shows the distribution of the 
full AGES parent sample (small black points) and the 
PRIMUS basis set. The size and color of the symbols rep- 
resent the basis weight of each template (redder/larger 
symbols show basis galaxies which characterize more par- 
ent sample galaxies at PRIMUS resolution) . In both the 
D4000 versus [O II] EW and color-magnitude parameter 
space, basis galaxies clearly dominate the blue/star form- 
ing portions of the diagrams; these galaxies show a large 
diversity in emission line strengths and continuum shapes 
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Fig. 15. — Random selection of 15 basis galaxies from the tem- 
plate used to fit PRIMUS spectra. Basis galaxies were selected by 
asking for the set of galaxies from fits to AGES objects which fully 
spanned the AGES galaxy sample if convolved to PRIMUS reso- 
lution. As the majority of the diversity in galaxy properties arises 
from emission line strengths and continuum shapes in blue galaxies, 
most of the PRIMUS basis galaxies are star-forming systems. 

compared to more quiescent galaxies. The red/quiescent 
basis galaxies, on the other hand, have very high basis 
weights; these basis objects describe a large fraction of 
the AGES parent set when convolved to PRIMUS reso- 
lution. 

5.2. AGN and Stellar Spectral Templates 

In addition to galaxies, PRIMUS observed broadline 
AGNs; Paper I includes the details of the AGN selection 
in PRIMUS, but in brief we often did nothing to remove 
AGN from our galaxy sample and in some fields we ac- 
tively targeted objects with color characteristic of AGNs. 
In order to measure redshifts for these objects and en- 



Fig. 16. — Properties of the PRIMUS basis sample (colored 
points) compared to the space spanned by the full AGES galaxy 
sample. In both panels, the color and size of data points shows 
the basis weight for each basis galaxy. Redder (larger) data points 
represent basis galaxies that can be used to represent more AGES 
sample galaxies when convolved to PRIMUS resolution. The vast 
majority of the variation in this basis set is in the blue/star forming 
galaxies. Left: Comparison of the strength of the 4000 A break, 
D4000, with the equivalent width of [O II]. Right: Coverage of the 
PRIMUS basis galaxies in restframe color-magnitude space. 



12 



sure that only galaxies are being assigned redshifts from 
galactic templates, we also fit each PRIMUS object with 
broadline AGN templates during redshift fitting. 

Obscured AGNs in which the optical spectrum has con- 
tributions from narrow emission lines from the AGN but 
is not dominated by broad-line emission are fit using the 
galaxy templates. The AGES parent sample contains 
examples of these Type II AGN and the stellar contin- 
uum dominates the broad shape of the spectra of these 
objects, and thus the galaxy templates are sufficient for 
determining redshifts. It should be noted, however, that 
due to the low resolution of PRIMUS, the presence of 
narrow emission lines from AGNs does not allow us to 
classify galaxies as AGNs. Throughout this section, we 
focus on creating a set of broad-line AGN templates used 
to classify objects which are dominated by the non-stellar 
AGN continuum and broad emission lines. 

As a basis for our set of AGN templates, we 
use the SPSS compos ite quasar spectrum from 
iVanden Berk et all (|2001l) . While this composite con- 
tains the broad average of quasar features, the full pop- 
ulation of AGNs have a wide range of spectral indices, 
a, and reddening due to dust. To span this range of pa- 
ra meters, we create a grid of AGN models based upon 
the IVanden B erk et aTT([2001l ) composite but for which 
we have modified the spectral index of the AGN contin- 
uum to values of a = -2.5,-2.0,-1.5,-1.0,-0.5,0.0 
and have included dust reddening with E(B — V) = 
0.0, 0.05, 0.1 assuming a dust model similar to the SMC. 
The range of E(B — V) values agrees well with the 
ob served range o f redde ning observed in SDSS quasars 
in iHopkins et al.l ([2004). We find that the population 
of AGN in the PRIMUS sample is well spanned by 
these modifications in a and E(B — V) and tests with 
broader ranges in each parameter showed the majority 
of PRIMUS AGN to be best fit inside the final range 
specified above. 

In addition to the variety of continuum shapes found in 
AGNs, as one observes AGNs at increasingly higher red- 
shifts, attenuation from the intergalactic medium (IGM) 
plays a critical role in shaping the rest-frame UV (ob- 
served optical) spectra. Neglecting this attenuation can 
lead to large differences between model spectra and ob- 
served AGNs at z > 3. With the depth probed by 
PRIMUS, we expect many quasars at these redshifts to 
be present in the sample, and thus we include a redshift- 
dependent attenuation from the IGM in our library of 
AGN template spectra based on an the redshif t evolu- 
tion o f the IGM attenuation curve presented in IMadaul 
(fl995l) . 

In addition to galaxies and AGN, the straightforward 
flux limit cuts in the PRIMUS target selection include 
galactic stellar sources as well, especially at faint mag- 
nitudes where photometric separation between stars and 
galaxies becom es more difficult. We utilize the empirical 
iPicklesI (|1998l ) stellar library to create a sample of stellar 
templates which we compare to each observed PRIMUS 
spectrum in order to determine if the source is stellar in 
nature. Due to the low resolution of PRIMUS spectra, 
subtle changes in the stellar absorption features or spec- 
tral shape are not important for our application; the goal 
with these stellar spectra is to ensure that stellar objects 
are not misclassified as galaxies when we fit redshifts. 




2000 4000 6000 8000 10000 

woveleogth (A) 



Fig. 17. — Example of a PRIMUS fit to an object identified as 
a 2 = 0.62 galaxy. The top panel shows the observed PRIMUS 
spectrum from both slit A (black) and slit B (blue). The dashed 
line shows the best-fitting model convolved to PRIMUS resolution. 
In the bottom panel, the grey line shows the full-resolution best- 
fitting galaxy model. The data points with errorbars mark the 
broadband UV and optical photometry available for the object. 

5.3. Fitting Method 

With the spectral template sets constructed, we per- 
form redshift fitting of each PRIMUS object. Through- 
out the fitting process, we fit the slit A and slit B observa- 
tions from our nod & shuffle method simultaneously. In 
addition to the PRIMUS spectra, we also simultaneously 
fit the observed optical photometry of each source. The 
addition of the photometry both adds information out- 
side the spectral window of the PRIMUS spectroscopy 
and allows us to find the best fit in the possible pres- 
ence of long- wavelength spectrophotometry errors. As 
we have not corrected for the large (and highly uncer- 
tain) slit losses in our PRIMUS spectra, we do not force 
the simultaneous fitting of the photometry and spectra 
to agree in absolute normalization. 

When fitting galaxies, we first fit each object on a 
coarse grid with Az = 0.01 spacing from z = to z = 1.2. 
AGN are fit from z = 0.0toz = 5.0 with 181 grid points 
such that the ith grid point, zagn.j, is 

l + ^AGN,i = 1.0P. (14) 

At each grid point, we calculate the x 2 between the 
measured photometry (with associated errors) and the 
PRIMUS spectra. When fitting, we do not give the pho- 
tometry more weight than is given by the associated 
errors; in essence, photometry is adding ^6 more data 



13 



points of information to the 300 data points available 
from spectroscopy and thus does not dominate the fit. 
Due to possibility of long- wavelength spectrophotometry 
errors in the PRIMUS spectra, we fit each object with a 
model of the form 

M (A) = a M (A) (a^(A) + a 2 £ 2 (A)) (15) 

where Mo(A) is the original template spectrum and £(A) 
is a refluxing vector that account for possible errors in 
the spectrophotometry. In detail, £(A), is a term that 
characterizes a linear term in logA 



1 + £(A) = 2 



log (A/A m in) 
l°g (Amax/A m in) 



(16) 



When performing the fits, we place priors of 10% on a% 
and a 2 ; while we expect some level of spectrophotomet- 
ric errors on large- wavelength scales, we do not want the 
fitting to be dominated by driving the refluxing terms to 
unreasonable levels to compensate for actual differences 
between the selected model and the observed spectrum. 
This refluxing only applies to the comparisons between 
the PRIMUS spectra and convolved model. The pho- 
tometry are not fitted with any form of refluxing. 

Based on these fits we calculate the \ 2 between each 
template spectrum and the observed photometry and 
spectroscopy. For stars, we choose the best stellar fit 
to be the stellar model with the lowest % 2 . For AGN, 
we choose the best fitting AGN template at each grid 
point by choosing the lowest \ 2 model. For galaxies, 
as we have information about the real frequency of each 
galaxy template in the full AGES parent sample. Rather 
than choosing the lowest \ 2 model, denoted Xmim at eacn 
grid point, we calculate the effective \ 2 ■> Xeff> based on 
the basis weights, hi 



Xeff 



Xmin 



2 log ^ bj exp 



x 2 

A. mm 



(17) 



For galaxies and AGN, we find the minimum in \ 2 
versus redshift in our grid. We then fit a finer linear grid 
with Sz — 0.001 in an interval of \z — z m { n \ < 0.03 around 
the minimum found in the course grid, z m ; n . Based on 
this finer grid, we calculate the final redshift of the galaxy 
or AGN as 



^bcst 



z exp 



xW 2 



x 2 - 

A mm 



dz 



(18) 



and the associated error is the second moment of the 
probability distribution 



z 2 exp 



Xn 



dz. 



(19) 



Once we have fit each PRIMUS object with galaxy, 
AGN, and stellar models, we must choose which of these 
models is the best representation of the measured data. 
If the final star or galaxy x 2 1S the lowest among the 
models, we choose to classify the PRIMUS object as a 
star or galaxy. After manually inspecting objects with 
Xaqn < Xgai' including objects with known redshifts and 
classifications from DEEP2 and zCOSMOS, we find that 



a threshold of Xagn + 100 < X ga i provides correct red- 
shift identification for objects that are clearly AGN (ei- 
ther from inspection of PRIMUS spectroscopy or previ- 
ous high- resolution spectroscopy). 

Figures QUEH and [H show example PRIMUS fits. In 
each plot, the top panel shows the PRIMUS spectrum ob- 
served through slit A (black) and slit B (blue) as well as 
the best-fit spectrum convolved to PRIMUS resolution 
(dashed line). The data bars on the PRIMUS spectra 
are counting statistics and do not include possible sys- 
tematics from sources such as flat fielding or errors in 
our scattered light model. The bottom panel shows the 
best-fitting high-resolution model as well as the broad- 
band photometric data points available for each object 
(data points with errorbars). 

6. REDSHIFT ACCURACY 
6.1. Quality Control 

Before using the PRIMUS redshifts for scientific work, 
it is vital to remove objects which could be contaminated 
due to systematics from the extraction or redshift fitting. 

The majority (90% see Paper I) of PRIMUS redshifts 
do not change when photometry is included in the fit, 
however the addition of photometric measurements to 
the PRIMUS fits does improve the spectroscopic fits of 
the remaining 10%, especially at fainter magnitudes. We 
thus require that any robust redshift measurement from 
PRIMUS include at least three bands of photometry. 
While some of the redshifts for objects with fewer than 
three bands of photometry may be correct, we find the 




Fig. 18. — Same as Figure [T7l but an example of a PRIMUS fit 
to an object identified as a star. 
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FlG. 19. — Same as Figure fTTl but an example of a PRIMUS fit 
to a z = 1.67 AGN. 
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FlG. 20. — Comparison of PRIMUS redshift with redshifts 
obtained from higher-resolution spectroscopy from COSMOS, 
DEEP2, and WDS for Q = 4 quality rating. The bottom panel 
shows the histogram of the difference between the PRIMUS and 
high-resolution redshift for each survey as well as the full compari- 
son sample. Overall, we find a ~ 0.5% dispersion between PRIMUS 
and the high-resolution sample with ~ 8% of the galaxies falling 
more than 5<r from agreement. 



failure rate for objects with less than three bands of pho- 
tometry is dramatically higher than for objects with more 
bands. Similarly, we flag any objects with more than 70% 
of the fitted wavelength range masked in the PRIMUS 
spectroscopy. These objects are predominately objects 
that fell near the edges of CCDs on the IMACS focal 
plane. Clearly, missing large sections of the object spec- 
trum limits our ability to measure a robust redshift for 
these objects. 

We further flag objects if the fitting of the spectrum 
shows signs of possible problems. If the best fitting red- 
shift was at the extremes of the redshift range fit z = 1.2 
for galaxies or z — 5.0 for AGN the object is flagged 
as a possible problem. If the best fitting redshift is at 
the extreme of the fitting range then we are unable to 
robustly measure a maximum in the probability distri- 
bution function. Additionally, we flag galaxies if the best 
fitting galaxy redshift is less than a 5cr (A% 2 < 25) de- 
tection compared to other redshift minima in the grid. 
Typically, these are low signal-to-noise objects with very 
flat probability distribution functions in redshift and the 
minimum is a noise spike. Lastly, we flag any objects 
with x 2 > 10 4 ; these are objects for which no model in 
our library performed well in fitting the object and thus 
no redshift information from these objects is trusted. 

The final test we perform checks the total available 
information available in the measured photometry and 
spectroscopy. We fit each object with a grid of pow- 
erlaw models with indices running from a — —3.0 to 
a = 3.0. If our spectra and photometry provide little 
information beyond a simple powerlaw model with no 



spectral features, we do not trust the measured redshift 
for the object. If the best fitting powerlaw has a x 2 °f 
xLt - Xpowcriaw > ~ 10 ( Le the best nttin S template 
model isn't better than ~ 3cr compared to a simple pow- 
erlaw) we flag the object as suspect and remove it from 
our final statistical sample. 

Once we have a clean sample of PRIMUS galaxies we 
assign each object a numerical redshift quality value. For 
each galaxy, we construct the statistic, £, which quanti- 
fies the separation of multiple minima in redshift space 
for each galaxy and AGN in the sample. Explicitly, we 
measure the difference in depths of the best and second 
minima in \ 2 space, A% 2 



Ax 2 = Xand 



Xbcst 



and calculate 



C = 1000 



(cr 2 /(l + Zbest)) 



(Ax 2 ) 



2\V2 



(20) 



(21) 



The numerator quantifies the width of the best peak in 
the P(z) distribution and the denominator quantifies the 
separation, in a statistical sense, between the goodness 
of fit in the first and second minima. In extremes, if an 
objects has a very broad peak in P(z) or if the difference 
in x 2 between the first and second peaks are very close, 
( increases. We define three confidence classes based on 
this statistic. 



f 2 if log C > 0.3 
Q = \ 3 if -0.3 < logC < 0.3 
{ 4 if log C < -3.0 



(22) 
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Fig. 21. — Same as Figure l20l but for galaxies with Q = 3 and 
Q = 4. The inclusion of Q = 3 adds some dispersion compared to 
the Q = 4 sample - the Q > 3 sample shows a lcr dispersion of 
0.7% and 15.3% rate of outliers more than 5a from agreement. 

The thresholds for each class are empirically determined 
by finding values which maximized sample completeness 
while minimizing outlier rate. With this classification, 
objects with Q = 4 have the most certain redshifts. For 
most statistical studies, a combination of Q = 3 and 
Q — 4 provides a sample of galaxies with confident red- 
shifts and low outlier rate. Throughout the rest of this 
work, we utilize a Q > 3 sample when measuring survey 
completeness. Table Q] lists the scatter and outlier rate 
for all three quality classes as well as the fraction of the 
PRIMUS primary galaxy sample which were classified 
into each bin. 

6.2. Redshift Success 

In addition to the science objects observed as part of 
PRIMUS, we also observed targets with known redshifts 
in DEEP2 and VVDS fields for calibration purposes. Af- 
ter our initial observations, zCOSMOS released an early 
redshift catalog which we further use to create a sample 
of PRIMUS objects with known redshifts from higher- 
resolution spectroscopic observations. In order to quan- 
tify the redshift success from our PRIMUS observations 
and redshift fitting, we compare our best-fit PRIMUS 
redshifts to this sample of known high-resolution red- 
shifts. 

Figure HOI shows the best fitting PRIMUS redshift com- 
pared to higher-resolution spectroscopic redshifts from 
DEEP2, zCOSMOS, and VVDS for Q = 4 objects and 
the associated distribution of differences. Overall we find 
excellent agreement between the two redshift measure- 
ments with a ~ 0.5% scatter in Az/(1 + z). Defining 
catastrophic outliers as objects with Az/ (l + z) > 5a we 
find a total outlier rate of 8%. Similarly, Figure |2"T1 shows 
the same distribution for all Q = 3 and Q = 4 galaxies. 
We find a 0.7% scatter rate and 15% outlier rate. The 



Q > 3 sample is the best suited for statistical studies of 
galaxy properties as it is more complete. 

One of the largest sources of error in PRIMUS red- 
shifts is the presence of degeneracies between two model 
spectra at different redshifts. At PRIMUS resolution, a 
galaxy at z > 0.7 with a strong 4000 A break can be dif- 
ficult to differentiate from a bluer galaxy with a Balmer 
break and a strong [O II] . emission line at slightly higher 
redshift. The addition of high-quality photometry can 
help break this degeneracy, especially in the U-band. As 
an example of the aid that U-band photometry can pro- 
vide to our fits, in regions with U-band photometry, we 
find a 25% tighter relation between high-resolution red- 
shift and PRIMUS redshift and the number of outliers 
located more than 5<r from the one-to-one relationship 
dropped by a factor of two. 

7. SURVEY COMPLETENESS 

Regardless of the effort invested into any wide-area, 
deep, galaxy redshift survey, 100% of the objects of in- 
terest will not yield redshifts. incompleteness can arise 
at many phases of the survey process; galaxies may not 
be present in photometric catalogs due to catalog depth, 
technical limitations or survey time limitations may pre- 
vent the targeting of every galaxy selected as a possible 
target from photometric imaging, and finally observed 
spectra may fail to provide a reliable redshift measure- 
ment. Proper use of data obtained from galaxy redshift 
surveys requires the careful correction for these sources 
of incompleteness in order to not bias statistical anal- 
ysis of the sample. As the photometric catalogs we 
target from in PRIMUS are often significantly deeper 
than the PRIMUS spectroscopic depth, we assume the 
incompleteness in the photometric catalogs to be neg- 
ligible and focus on targeting and redshift incomplete- 
ness in the PRIMUS observations. When selecting tar- 
gets for PRIMUS spectrosco py, we app li ed tw o a-priori 
sparse sampling criteria (see ICoil et al.l (|2011[) for more 
details). In order to ensure that the faintest galaxies, 
which dominate the photometric catalog, do not domi- 
nate the PRIMUS target selection, we define two bins in 
flux. Galaxies more than 0.5 mag brighter than the flux 
limit of the PRIMUS selection were sampled at a rate 
of 100% while galaxies below this limit were sampled 
at a rate of 30%. As this sparse sampling is performed 
during target selection and the sparse sampling weights 
were saved, we can correct for the magnitude-dependent 
sparse sampling exactly. Secondly, PRIMUS target se- 
lection utilized a density-dependent sparse sampling ap- 
proach to ensure that galaxies in (projected) dense re- 
gions were not undersampled. We can upweight each 
galaxy by the inverse of the sparse-sampling weight to 
reconstruct the full magnitude-limited sample. 

When designing slit masks for PRIMUS observations, 
not every object selected as a PRIMUS target can be 
observed due to slit collisions between PRIMUS targets 
themselves and between targets and higher priority tar- 
gets such as alignment stars and flux calibration stars. 
The effect of these slit collisions is minimized by observ- 
ing two slit masks at each PRIMUS field center and by 
using density-dependent sampling, but this does not al- 
low us to observe every possible PRIMUS target. In or- 
der to correct for missing targets from slit collisions, we 
calculate the number of PRIMUS galaxies which were 
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Fig. 22. — Redshift success as a function of color and magni- 
tude for each of the PRIMUS fields. Contours enclose 25%, 50%, 
75%, 90%, and 99% in the targeted PRIMUS sample while the 
colored background shows the fraction of targeted galaxies with a 
well-measured (Q > 3) redshift in PRIMUS. Due to the inhomo- 
geneous filter systems utilized in the optical imaging across each 
field, we utilize synthesized r — i color and i-band magnitude based 
on kcorrect fits to the observed photometry in this comparison. 
Overall, the primary factor in our redshift success rate is the i-band 
flux of the object; little trend is seen with respect to an object's 
color. The average redshift success is > 75% brighter than i = 21 
and declines to ~ 45% at i ~ 22.5. 

observed and divide by the number of PRIMUS galax- 
ies available as a function of the number of masks each 
object had the opportunity to be observed on and the 
number of potential conflicting galaxies, N c . Averag- 
ing over all of the PRIMUS fields except XMM-SXDS, 
our completeness in areas observed by two masks is 97% 
for Nc < 2 and declines with increases Nc to 82% for 
N c = 6. In areas covered by four masks (in regions of 
the COSMOS field), our average completeness is 98% for 
1 < N c < 4 and declines slightly to 94% for N c = 6. In 
XMM-SXDS, we observed a factor of two more F-stars 
per mask than our other fields, which leads to lower over- 
all completeness (96% for N c = 2 and 75% for N c = 6). 

Finally, we must correct for galaxies which were ob- 
served by PRIMUS but for which we were unable to mea- 
sure a reliable (defined as Q > 3, here) redshift. In order 
to quantify this incompleteness and explore its depen- 
dence on galaxy type, we count the number of observed 
galaxies as a function of both observed magnitude and 
color. As the photometric systems between each of the 
imaging datasets used for target selection used different 
filter sets, the exact filters used for this comparison are 
different between fields. In general, we choose either the 
i or R bands when quantifying trends with magnitude 
and choose a color with a large lever arm across the op- 
tical spectrum (typically g — i or B — R) when looking at 
trends with observed-frame color. We then derive the fi- 
nal spectroscopic weight in this binned color-magnitude 
space. We define the final weight as the total number 



of galaxies in a bin with robust redshift measurements 
divided by the effective number of galaxies (that is the 
total of the PRIMUS statistical weights) from the full 
primary sample in the bin. We use this method to en- 
sure that galaxies from populations with few successful 
redshift measurements are still given full weight in the 
final primary sample. Figure [5U shows the complete- 
ness for each of the PRIMUS fields. In each panel the 
contours show the concentration of objects in the color- 
magnitude plane and denote 25%, 50%, 75%, 90%, and 
99% of the population. The color denotes the fraction 
of observed PRIMUS objects which received a reliable 
redshift measurement. In general, the flux of the object 
is the dominant variable when describing the PRIMUS 
redshift incompleteness - little variation is seen in the 
completeness as a function of color at fixed flux. The 
average redshift success is > 75% for galaxies brighter 
than i = 21 and declines to ~ 45% at i = 22.5. 

8. CONCLUSIONS 

PRIMUS has completed observations of 185,105 galax- 
ies with a low-dispersion prism over 9.1 deg 2 suitable for 
statistical investigates of galaxy evolution since z *~ 1. 
In our primary sample of galaxies with redshift mea- 
surements suitable for statistical work, PRIMUS in- 
cludes 130,000 galaxies with typical redshifts errors of 
er 2 /(l + z) — 0.005 when compared to higher-resolution 
spectroscopic redshift measurements. PRIMUS is the 
largest faint galaxy survey completed to date. The high 
targeting fraction (~ 80%) and large survey volume allow 
for precise measur ements of galaxy properties and large- 
scale to z~ 1. In iCoil et al.l (| 2 1 If ) . we summarize the 
motivation, observational techniques, target selection, 
slitmask design, and provided details of the PRIMUS ob- 
servations. In this paper, we detail the data processing, 
spectral extraction, redshift fitting, and survey complete- 
ness of the PRIMUS. 

The dramatic increase in multiplexing ability added 
when using the low-dispersion prism technique developed 
for PRIMUS compared to more traditional spectroscopic 
techniques at higher-resolution allows for a very efficient 
mapping of large areas of the sky using only a fraction 
of the observing time. The value of high-efficiency spec- 
troscopic follow up will increase as the next generation 
of large area deep sky surveys such as the Dark Energy 
Survey and LSST begin providing public data. The vol- 
umes probed by these imaging surveys will be impossible 
to fully spectroscopically probe using traditional tech- 
niques. 

While the prism technique has many strengths, espe- 
cially in terms of survey efficiency, it has several limi- 
tations; the effect of these limitations can be minimized 
through proper planning and careful survey design. As 
much of the spectral information on small scales (ab- 
sorption and emission lines features) will be unresolved 
with prism spectroscopy, the more information available 
to break degeneracies in spectral fitting will result in 
more precise redshift measurements. For example, in 
PRIMUS, a spectral break at z > 0.7 may arise due to a 
strong 4000 A break in a red galaxy or a Balmer break 
with an unresolved [O II] emission line in a star forming 
galaxy. The addition of external photometry, especially 
in the [/-band, helps to breaks this degeneracy. Having 
broadband photometry which spans the optical spectrum 
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provides significant improvement to redshift fits. 

Traditional higher-resolution spectroscopic surveys of- 
ten utilize sky fibers, extended slits, or nod & shuffle pro- 
cedures to characterize the sky emission in target spectra 
as the signal from targets of interest are many orders of 
magnitude fainter than the sky and thus any contami- 
nating emission not removed will bias the final extracted 
spectra. This problem is magnified with low dispersion 
spectroscopy as the entire trace of objects of interest may 
only be a few hundred pixels of information. Small scale 
errors that may affect only a small portion of the spec- 
trum in high-resolution spectroscopic work (and which 
can easily be masked out when fitting for rcdshifts) re- 
moves a significant fraction of the spectral information 
at lower-resolution. Robust sky subtraction is essential, 
though the most effective technique is strongly depen- 
dent on the instrument, goals, and depth of each survey. 

A key problem that required significant effort to cor- 
rect in our PRIMUS observations was contamination of 
our spectra from large-scale scattered light in the final 
images. This scattered light is likely insignificant for 
the majority of users of the instrument as it exists on 
large enough scales to be removed when subtracting sky 
emission in direct imaging. The scattered light being de- 
pendent on the flux from nearby pixels means that most 
traditional high-resolution spectroscopy applications will 
likely be minimally effected as the counts per pixel are 
often significantly smaller than when working with low- 
dispersion spectra (as the sky emission is resolved rather 
than being unresolved in the red). This underlines the 
critical importance in understanding the instrument and 
optical system when performing low-resolution spectro- 
scopic work. Often corrections unneeded, or even unno- 
ticed, in more traditional applications will play an im- 
portant role in the final success of the survey. 



Once a characterization of the available photometric 
data, a coherent plan to remove contaminating emission 
from object spectra, and a firm understanding of the in- 
strument and optical systems involved have been estab- 
lished, low-resolution spectroscopic surveys provide very 
efficient means to observed thousands of objects simul- 
taneously. Early success with PRIMUS led to a number 
of current prism surveys to probe different areas of pa- 
rameter space including the high redshift universe, the 
nature of galaxy cluster members, and the properties of 
colliding galaxy clusters. 
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TABLE 1 

PRIMUS Redshift Confidence Classes 



Class 


°\5z/(l+z) 


Outliers a 


Sample Fraction 13 


4 


0.005 


7.85 


49.2 


3 


0.022 


5.32 


21.6 


2 


0.050 


5.06 


29.2 



a Fraction of objects with known redshifts deviating 
more than 5<r from agreement. 

b Fraction of PRIMUS primary galaxies which re- 
ceived the specified class designation. 



